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SYMBOLS 


j imaginary operator, 

s Laplace variable 

Sf Laplace transform operator 
t time, sec 

T computer step size, **frame time,” sec 
u{t) unit step at time t 
X continuous input variable 

X Laplace transform of x 

x^ sampled input variable 

X* Laplace transform of x* 

y piecewise continuous DAG output variable 

Y Laplace transform of y 

e error between x and y 

T continuous time variable in interval T 

0) frequency, rad/sec 

(*) time derivative 

definition 
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FIRST-ORDER-HOLD INTERPOLATION DIGITAL- TO- ANALOG CONVERTER 


WITH APPLICATION TO AIRCRAFT SIMULATION 
William B. Cleveland 
Ames Research Center 


SUMMARY 


Those who design piloted aircraft simulations must contend with the 
finite size and speed of the available digital computer and the requirement 
for simulation reality. With a fixed computational plant, the more complex 
the model, the more computing cycle time is required. While increasing the 
cycle time may not degrade the fidelity of the simulated aircraft dynamics, 
the larger steps in the pilot cue feedback variables (such as the visual scene 
cues), may be disconcerting to the pilot. 

The First-Order-Hold Interpolation (FOHI) Digital-to-Analog Converter 
(DAC) is presented as a device which offers smooth output, regardless of cycle 
time — a significant improvement over the conventional Zero-Order-Hold (ZOH) 
DAC and the First-Order-Hold Extrapolation (FOHE) DAC. The Laplace transforms 
of these three conversion types are developed and their frequency response 
characteristics and output smoothness are compared. The FOHI DAC exhibits a 
pure one-cycle delay. Whenever the FOHI DAC input comes from a second-order 
(or higher) system, a simple computer software technique can be used to 
compensate for the DAC phase lag. When so compensated, the FOHI DAC has 
(1) an output signal that is very smooth, (2) a flat frequency response in 
frequency ranges of interest, and (3) no phase error. When the input comes 
from a first-order system, software compensation may cause the FOHI DAC to 
perform as an FOHE DAC, which, although its output is not as smooth as that of 
the FOHI DAC, has a smoother output than that of the ZOH DAC. 

INTRODUCTION 


Over the last decade, the field of real-time, man- in- the- loop aircraft 
simulation has undergone some basic changes. The parallel computing ability 
of analog computers has given way to the serial digital computer with its 
obvious computational advantages and some not so obvious disadvantages. Early 
digital computers had insufficient speed to simulate an aircraft, control 
cockpit simulators, and acquire data (see fig. 1), all without producing a 
**ratchetlike’* response in the simulation. Now, however, digital computers are 
so fast that the serial computation of the kinematic equations in real-time is 
no longer a serious problem, but discreteness in the outputs remains a problem 
for simulation. From a psychological standpoint, the discreteness of the data 
output from the computer can be distracting to the pilot attempting to fly a 
simulated aircraft. He expects to see smooth continuous changes in his flight 
information, just as he would in the real aircraft, but the discreteness of 



the computer output may degrade the 
realism of the simulator to the 
point where he "flies the simulator" 
rather than the aircraft. Experi- 
enced pilots can accommodate to or 
compensate for many of the indica- 
tions that they are flying a simu- 
lated aircraft rather than a real 
aircraft ; however, it is distract- 
ing to the pilot to contend with 
spurious flying cues that are due 
to output data discreteness and not 
to the simulated aircraft 
characteristics . 

There are two causes of the 
discreteness of data output from 
These are quantization effects 
that are due to the number of bits in the DAC word and the effects that are 
due to the discrete nature of the output. Figure 2 illustrates the latter. 

The roughness due to DAC quantization is alleviated by increasing the number 
of bits in the DAC word length. However, the roughness due to time discrete- 
ness is a complex phenomenon which involves computer hardware and software. 

The problem for the simulator pilot is most acute when the instrument dynamics 
can closely follow the stairstep-type signal or when the steps are magnified 
in some way, such as with angular displacements in a visual scene. A typical 
visual scene is modeled to scale and scanned by a television camera on a 
servo-driven gantry which responds to an angular pitch "up" command by shift- 
ing the scene downward on the television monitor. The angular discreteness is 
translated into a positional discreteness on the visual scene monitor through 
a multiplication of the distance from pilot's eye to the monitor. Reduction of 
these jumps is desirable to avoid pilot distraction in a research environment. 

Many roughness problems can be cured by increasing computation speed, but 
most complex real-time programs are already running near the maximum speed of 
the computer. This conflicts directly with the need to lengthen the cycle 
time to allow more complex (and therefore more realistic) models to be simu- 
lated. The FOHI DAC, analyzed within this report, copes with this problem by 
eliminating or reducing the discreteness in DAC outputs to the simulator. 

DAC's are a rather small component in a total piloted simulation system 
which includes a computer, a cockpit, and simulators. However, they are the 
important link within the pilot's closed-loop control system. In the analysis 
of any control loop, the model of each element is needed to assess the whole. 
Therefore, the frequency response and the discreteness effects of two well- 
known types of DAC's, the ZOH and FOHE DAC's, are compared with those of the 
FOHI DAC in the following discussion. 



Figure 1.- Elements of the piloted 
digital aircraft simulation. 


the digital computer via the ordinary DAC. 
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x(t) 


(a) Continuous function 





Figure 2,- Comparison of various DAC outputs with a corresponding continuous 

signal . 


DIGITAL- TO-ANALOG-CONVERTER FREQUENCY RESPONSE 


The function of the DAC is to change the form of a variable from a pat- 
tern of bits in the digital word into a continuous (normally piecewise con- 
tinuous) analog voltage signal. Figure 2 compares the ZOH, FOHE, and FOHI DAC 
output signals versus a continuous counterpart signal. By far the most common 
type in use today is the ZOH DAC. Figure 2(b) shows the way a value is held 
constant between sample points which occur at interval T, 

In an effort to obtain a smoother output from the digital computer, 
First-Order Hold (FOH) DAC * s may be employed. Figure 2(c) illustrates the time 
history of an FOHE DAC, which uses one past in addition to the present value 
of discrete input data to extrapolate a more representative output than that 
of the ZOH DAC. Inspection of figures 2(b) and (c) shows that, as the fre- 
quency content of the continuous signal increases and as T is increased, the 


3 



discontinuous jump in the outputs also increases. If interpolation, instead 
of extrapolation, between sample points is used on the FOR DAC, there are no 
discontinuous jumps, since the linear portions of the output start and end on 
sampled values, as seen in figure 2(d). 

The transfer function of a DAC is derived next in order to make meaning- 
ful comparisons between the three types of DACs discussed. 


General Transfer Function of Linear Hold DACs 

For zero- and first-order DAC’s, the output y{t) is a piecewise linear 
function of the input and 


x*(t) 

x(t) ^ DAC y(t) 

T 

The output function is a constant or ramp in the nth T period and can be 
represented by 


DAC 


y„(x) = A^r + 0 < T < T (1) 

SO that 

00 

y(t) = ^ ( a^t + s„) (2) 

n=-oo 

where 

^ ~ 1)^] + (3) 

5^ = - 1)T] 4- (4) 

and where and are constants determined by the particular DAC. For 
example, the ZOH output has a constant value equal to x{nT) throughout the 
interval t. Thus a-^ - ^ ^2 “ 1.0. 

Taking the Laplace transforms, one obtains 


“ i«(n+l)T 


OO M 

I(s) = Sf[yit)] = ^ J 

n=-oo nT 


(A^t + dt 


(5) 


By letting t - t - nT, the above equation becomes 

^(3) - 1 / 

n=-oo \ •'(9 

- E 


T T 

-snT r -ST 7 . □ -snT “ 

xe dx + By^e 


f -ST 

j " 


6?T 


n=- 


4 




where 


5, = - (1 - - sTe~^^) 

* o 

and 

^2 = 7 

By using the shifting properties (discussed in the appendix), for 5^, one gets 

^ ia^x[(n - 1)T] + a^orCnT) 

n=-oo >^=-00 

mm G’V 

= X*(s)S'2(aje + 02 ) 

Similarly, for ^ 2 , we have 

^ = X*(s)S2(.b^e~^^ + b^) 

n=-oo 


The desired transfer function then becomes 

y(s) 1 , -sT . W-, ^ -sT^ . 1 -sT . , X -sTx 

{ci^e + a^) (1 - e ^ ^ 7 ■*■ ^2^ (1 “ ^ ) (6) 

The specific cases of the ZOH, FOHE, and FOHI DAC’s may be obtained by 
proper selections of the a and h constants. 


Zero-Order-Hold DAC 


The time function y (t) is seen, from sketch (a), to be a constant over 
the period T. In which case. 


y(t)=o:(nT) for nT < t < (n + 1)T . (7) 

y(t) 

Thus, comparing equation (7) with (1), * 

(3), and (4) it follows that 1 

x(nT) 1 

and Z ?2 “ 1 • I 

By substituting the derived values ' ' ' 

- , r , , - nT (n+1)T 

for and into the general trans- 
fer function, we have Sketch (a) 


Y(s) ^ 1 - 
Y*(s) s 


Jt 


( 8 ) 
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As shown in the appendix, given X(s) is the Laplace transform of x(t) y then 
the transform of the sampled variable x*(t) is 




-w'E x(s-J 2|I) 

n=‘“Oo > ' 


( 9 ) 


n=-o° 

The transfer function from x to y for the ZOH DAC may then be written 

-sT °° 


^(s) 

X{s) 


1 - e 


sT 


E ^ 


While the transfer function from a: to y is complicated by the repeating spec- 
trum, the multiplier 1/T renders the transfer function nondimensional and 
allows determination of the magnitude of the frequency response, as a function 
of the nondimensional product o)T. 


Thus , 


i r 1 ^(g) = 1 - g 

T z'- ’ T X*(s) sT 


-Ts 


( 10 ) 


By substituting s = joj, one may obtain the magnitude and phase of the sinu- 
soidal frequency response, as follows: 




uiT 


sin 




^ - - 




( 11 ) 

( 12 ) 


First-Order-Hold Extrapolation DAC 
From sketch (b) , one can see that 

^(t) = |a:(n2’) - a;[(n- l)r]| ^+x(nT) 

(13) 

for nT < t < (n + 1)T and 0 < t < T. 

Comparison of equation (13) with (1) , 

(3), and (4) yields 


a. 


= -«! = y ; = 0 and = 1 



which may be substituted into the general transfer function, equation (6) , to 
obtain 
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( 14 ) 


For the frequency response characteristics, the magnitude relationship is 


and the phase is 


T 


4(1 + (o2t2)1/2 . , 0)7 

= sxn^ — 


aj2r2 


^ ~ ^ oiT - (jiT 


First-Order-Hold Interpolation DAC 

Sketch (c) illustrates the time shift between input and output, which 
yields the time function 


y(t) = {x(nT) - x[(n - l)r} ^ 

+ cc[(n - 1)T] (17) 

for nT < t < (n + 1)T and 0 < t < ^. 

Comparison of equation (17) with (1) , 
(3), and (4) yields the constants 


x[(n-1)T] 


y(t) 


x(nT} 



^ 1 
T 


= -a 1 = &^ = 1 and = 0 


nT (n+1)T 

Sketch (c) 


which, upon substitution into the general transfer function, equation (6) , 
yields 


x*{s) 




Since 


f <=FI<°> ■ [f 

the magnitude and phase relationships are 


IJ1 (t? ^ ) 


4 . 2 

2m2 2 


o)^r 


Gp^Ux) = -(or 


(15) 

(16) 


(18) 

(19) 

( 20 ) 
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DAC SMOOTHNESS 


A measure of smoothness is the maximum error between the DAC output and a 
corresponding continuous input sine wave signal x = A sin h)t. For the ZOH, 
the error is 


e = A 



- A s±n (nwT) 


0 < T < T 


( 21 ) 


which is a maximum near the zero crossing of the sine wave. If we let the 
time nT be at the crossing, the maximum error occurs when t = T. The maximum 
error is then given by 


^max = A s±n oiT 


and, for small coT, the maximum error is 


^max/^ ‘^2’ (22) 

The maximum error for FOHE occurs when cot is near tt/2. The FOHE error 
relation is 

e = ^sinj^l^^ ~ ^ sin(nco^?^) + ^ ^ sin[(n- l)coT]| 0< t <7 

(23) 

If the time nT occurs at the peak of the sine wave, then n = Tr/2oo!r. The maxi- 
mum error occurs at t = (n + 1)!T. On substituting these values in the FOHE 
error relation, the solution for the maximum error is given by 

l^maxi/A (24) 

The maximum error for FOHI also occurs when u)T is near tt/2. Assuming the 
sample points are equally spaced about the peak amplitude, the maximum error 
occurs at the peak of the sine wave. Thus, the FOHI error relation is 

e = y - |i4g^j^(no)r) + A ^ sin[(n + 1)0)7] - A — sin(no)7)j 0 < t < 7 

(25) 

Under these conditions, nT occurs at t - -n/lu} - T/2^ so that n = TilloyT - 1/2, 
The maximum error occurs when x = 1/2. On substituting these values in the 
FOHI error relation and assuming small o)T, the solution for maximum error 
provides 


^max _ 0)^^^ 


(26) 
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The results are as one might expect from inspection of figure 2; the ZOH DAC 
exhibits the most error and the FOHI DAC the least. Table 1 gives a compari- 
son of the relative DAC error for the ooT of interest in aircraft simulations. 


TABLE 1.- RELATIVE SMOOTHNESS USING THE MAXIMUM ERROR IN THE DAC 
REPRESENTATION OF A SINE WAVE OF AMPLITUDE A AND FREQUENCY cj 


uiT 

Relative error 

^max/A 

ZOH 

FOHE 

FOHI 

0.1 

0.1 

0.01 

0.00125 

.2 

.2 

.04 

.005 

.3 

.3 

.09 

.01125 

.4 

.4 

.16 

.02 

.5 

.5 

.25 

.03125 

.6 

. 6 

.36 

.045 


COMPARISON OF THE DAC CONFIGURATIONS 


Figure 3 illustrates the magni- 
tude and phase for each of the three 
DACs under consideration. The ZOH 
DAC appears to be superior in terms 
of less phase lag over the range 
shown and also in the way its magni- 
tude is relatively flat at low 
values of ooT. In some control sys- 
tems, (joT might range close to the 
maximum allowable of cjT = tt, whereas, 
in simulations of aircraft, such 
things as numerical instabilities of 
integration algorithms require 
0 )T < 0.6 in practical cases . For 
example, the second-order Adams- 
Bashforth predictor, an integration 
algorithm commonly used in simula- 
tion, is numerically unstable when 
ix)T = 1.0; the third-order Adams- 
Bashforth is unstable when cdT = 0.5. 

In this case, co is the natural fre- 
quency of the system under consider- 
ation. To obtain reasonable accuracy, 
uiT must be considerably smaller than 
the stability bound. 



Figure 3.- Magnitude and phase charac- 
teristics of zero- and first-order 
DAC ’ s . 


A second rationale for limiting coY to 0.6 concerns the aircraft itself. 

The value of co in aircraft simulations depends on the aircraft type: a trans- 

port, such as the DC-8, has its highest short-period frequency at approximately 
3.6 rad/sec; a high-performance aircraft, such as the A7, has its highest 
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Figure 4.- Enlargement of magnitude 
and phase characteristics of the 
DAC*s at low o)T. 


value at about 8.5 rad/sec; while a 
helicopter’s flapping frequency, such 
as the CH-46, has a value of about 
14 rad/sec. Using a T of 0.05, a 
representative value in aircraft 
simulations, is 0.7 for the CH-46, 
but for the fixed wing aircraft, wT 
is well below 0.6. Consequently, 
helicopter simulations present com- 
putation problems for digital com- 
puters and require special modelling 
techniques . 

j Figure 4 shows the phase and 

* magnitudes of the DAC’s over the 
lower range of toT. It is interest- 
ing to note that the FOHI and ZOH 
DAC’s are almost flat in magnitude 
over this range. The phase of the 
FOHE DAC is very good over this 
range, but near an o)T of 0.5, the 
DAC shows significant gain. 

The phase of the FOHI DAC is 
that of a pure one-T lag. If, 
somehow, a pure one-T lead could 
be introduced, the combination of 
lead and FOHI DAC would : be very 

smooth, have flat frequency 
g response over gjY' ranges of interest, 
and have zero phase error. 


Phase Compensation of the 

To show that the phase lag of the FOHI DAC 
consider a unit which might produce such a lag: 


FOHI DAC 

is equivalent to a one-i" lag. 


u(t) 



v(nT) 


Between V and u is the i’-lag device. In terms of Z-transforms , 

V{z)/U{z) = 


(27) 
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Now 


Z = = cos coT + j sin wT (28) 

so that upon substitution of equation (28) into (27) the magnitude and phase 
of the lag device are 


F(jg)) 

f/(jw) 


1 


(29) 


and 


7(ja))/t/(jo)) = -i^T (30) 

Thus, comparing equation (20) with (30) is seen that the FOHI DAG has the phase 
of a one-T time-lag device, A device which produces a pure one-T lead would 
be useful in compensating for the time lag, but, unfortunately, is not avail- 
able. However, one can resort to a computer programming procedure to intro- 
duce the desired one-T lead in certain applications. 

Consider the case in which a second-order system precedes the DAC. For 
example, in the visual display for an aircraft simulation, the camera servo is 
positioned in the translational direction X through a double integration of 
the X component of the aircraft ^s acceleration: 



On the digital computer, the integrations are performed numerically with an 
integration algorithm. As an example, the Euler integration method illustrates 
how a one-T lead may be generated through the double integration process. 

Applying the Euler algorithm in the double integration as illustrated 
above , we have 


x{nT) = i[(n - 1)T] + Tx[(n - 1)T] (31) 

x(nT) = x[(n - 1)T] + Tx[(n - 1)T] (32) 

The computational technique would normally be to solve equation (32) first 
for x(nT) and then equation (31) for x(nT ) , in that order. The DAC output at 
the time nT is x(nT) , Now, if the order of integration is reversed, the first 
calculation yields x(nT); the second integration is performed by using x(nT) 
to obtain :c[(n + 1)7]. When o:[(n + 1)7] is transmitted through the DAC at 
time n7, a full one-cycle lead has been obtained. When this technique is com- 
bined with the FOHI DAC, the phase lead of the computation cancels the phase 

lag of the DAC. 

Where there is no second-order (or higher) system preceding the DAC, the 
application of the FOHI DAC is less advantageous. However, in cases where 
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there is one integration preceding the DAC, it is an easy inatter to add a 
simple lead term, thus converting the FOHI to a FOHE DAC as shown in figure 5. 
The input to the DAC is of the form x + Tx^ where Tx is the lead compensation. 


1 + Ts 


— ~ (1 — 


1+Ts (i_e-Ts)2 
Ts2 


Ts2 


Lead 

+ 

FOHI DAC 


FOHE DAC 


Figure 5.- Relationship between FOHI and FOHE DAC’s. 
Phase Compensation for the ZOH DAC 


Since the ZOH DAC is widely used, it is worthwhile to describe one method 
of compensation for the phase characteristics of this DAC. If the first time 
derivative of the output is available within the digital computer, the addi- 
tion of fex to the normal output x will result in phase lead, which may be used 
to offset the phase lag of the DAC. In block form, the system is 


X*(s) 


1 +ks 



1-e-Ts 



s 



Y(s) 


Compensation 


DAC 


The input/output transfer function is 


7(s) _ 1 + ks _ -Ts . 
X*is) ~ s U e ; 


(33) 


The phase of Y(s) /X'^(s) is the phase of the compensation plus that of G^(s) 

, . - 1 7 cpy 

(j) = tan KW ^ 


When (j) = 0, 


For 0)7 < 0.6 


so that 


7 1 ^ 0)7 

k = — tan — 
0 ) 2 



k 


7 

2 
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While setting k - Til results in nearly zero phase, there is a minimal gain 
magnification with a value of 1.04 at an o)!? of 0.6. 


CONCLUSION 


The FOHI DAC offers a method of replacing the stairstep output of the 
common ZOH DAC with a smooth signal. Where a DAC is associated with a second- 
order (or higher) system, the programming procedure described will result in 
an output which is smooth and has a flat frequency response without phase 
error. By itself, the smooth signal is useful for cases where signal phase 
angle is unimportant, as found in many nonreal-time applications. In the 
real-time case, the DAC whose output commands a position servo system is often 
associated with a second-order system simulated in the computer. In the field 
of aircraft simulation, the visual simulators, in which a servo positions a 
camera over a model runway scene, is such a system. Conventional ZOH DACs 
produce ”ratchetlike” jumps in the scene when either the signal frequency con- 
tent or the computer stepsize is large. Through the use of FOHI DAC’s, the 
distracting ratchet response is eliminated, even for large step sizes. This 
technique effectively conserves computer power by eliminating the need to 
shorten computer step size to achieve needed smoothness. The time saved is 
therefore available to expand the simulation model for increased fidelity. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, CA 94035, May 27, 1976 
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APPENDIX 


TRANSFER FUNCTION OF A SAMPLED SIGNAL 


The sampled signal x*{t) of a continuous signal x(t) can be represented 
by a series of spikes (delta functions) multiplied against x{t) 


00 

= x{t) - kT) 


where 


( 1 t == kT 
0 else 

The Fourier Series representation of this summation is 


y) 6(t - kT) 
- 


oo 

^ exp(j27mt/T) 


where 


T/2 92^ . 

^ T I 6(t - nT)&x-p(-j2'mt/T)dt = 

-T/2 n=-°° 


so that the Laplace transform is 




M °° OO 

/ I S 

.OO y^xz—o 


exp {Qli\ntlT) 


-sT 

e 


dt 


oo 

IE ; 

^=— oo a 


3:(t)exp[ (-S + 


ImT 


Now if 


then 


x{s) 



/ jL \ “ Siy 

x{t)e 


dt 


x{s “ jlimit) = r x(t)e t727rn/T) 

•L.00 
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and so 


oo 

J^[x*(t)] = 1/T x(s - j2-im/T) 

n=-“ 


Alternative Version of the Transfer Function 

Whereas the first version of the transfer function yields a result with 
a repeating spectrum, this alternative is of more use in the mathematical 
operations dealing with the discrete sampled signals. Again, the sampled 
signal is represented by a summation of the sampled signals, but the "sifting” 
properties of the delta function are utilized rather than the spectral proper- 
ties of the Fourier Series approach. 




o 


x*(t)e dt 


oo ^oo 


/ 


x(t)6(t - kT)e dt 


= x{kT)e 


-skT 


It is convenient in working with discrete signals to use the following 
shifting property of the transform: 


OO 

£e{x[{k + n)T]} = ^ x[{k + n)T]e 


-skT 


oo 

= x(mT)e 


-smT 


^=-00 


nsT ^ . 

= 0 x*(s) 


NASA- Langley, 1976 
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